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Abstract We compute renormalization group fixed points and their spec- 
trum in an ultralocal approximation. We study of two competing 
non-trivial fixed points for a tliree-dimensional real iV-component field: the 
0(A^)-invariant fixed point vs. the cubic- invariant fixed point. We compute 
the critical value Nc of the cubic 0^-perturbation at the 0(A^)-fixed point. 
The 0{N) fixed point is stable under a cubic 0^-perturbation below Nc, above 
Nc it is unstable. The critical value comes out as 2.219435 < Nc < 2.219436 
in the ultralocal approximation. We also compute the critical value of at 
the cubic invariant fixed point. Within the accuracy of our computations, 
the two values coincide. 



1 Introduction 



Spin systems with an A^-component real field variable, governed by a Landau- 
Ginzburg type Hamiltonian, are of central importance both in Euclidean 
quantum field theory and classical statistical mechanics. The basic model 
is one with a global 0{N) symmetry, the invariance under a simultaneous 
rotation of all the spins. A prototypical Hamiltonial for such a model, say 
on a lattice A C Z^, is 



xeA I a=l \fi=l 



N 



It is important to study the influence of perturbations which explicitely break 
this symmetry to a smaller subgroup. For instance, in cubical crystals, one 
expects the spin interaction to react to the lattice structure. This suggests 
additional terms in the Hamiltonian which are not rotation symmetric, but 
invariant under the cubic group. The cubic group is composed of the per- 
mutations and reflections of the A^ components of the field. Such a cubic 
invariant term is 

N 

(2) 

xeA a=l 

When (|l|) is augmented by (0), a competition of renormalization group fixed 
points sets in, about which of them determines the long distance behaviour 
of the model. 

According to folklore, the 0{N) fixed point becomes unstable above a 
certain threshold value A^^. When this happens, the cubical fixed point dom- 
inates the long distance behaviour. The value of Nc is still being debated. 
Work on the e-expansion suggests that Nc > 3. This has been challenged 
by the studies |]12|, 0, ||, |[, which suggest that < 3. Recent Monte 
Carlo work |^ by Caselle and Hasenbusch again indicates that A'c should be 
very close to 3. Within their precision, this result is compatible with that of 
Kleinert and collaborators 0, |§, |^ . 

To our knowledge, the problem has not yet been looked at in the frame- 
work of the block spin renormalization group ||16[. To study the domains 



of attractions in a model with several competing fixed is a fundamental and 
challenging problem as it underlies Wilson's explanation of universality. A 
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notorious trouble in this business is that these Hamiltonians in general de- 
pend on an infinite number of interaction terms. To control fiows with a large 
number of couplings is a very difficult task. Furthermore, when symmetries 
are reduced, the number of couplings tends to proliferate dramatically. 

This paper contains a study of both the 0(iV)-invariant fixed point and 
the cubic fixed point in the framework of the hierarchical or ultralocal approx- 
imation to Wilson's renormalization group. Here we restrict our attention 
to the case of D = 3 dimensions. In the hierarchical approximation, the 
effective Hamiltonians are restricted to a (non-standard) kinetic term plus 
local interactions. The local interactions are given by a potential, which is a 
function of variables. 



2 Hierarchical Renormalization Group 

The hierarchical renormalization group is a block spin renormalization group 
for so called hierarchical spin models. They are spin models with a non- 
translation invariant kinetic term designed to make the renormalization group 
local. The hierarchical model and the full model belong to different universal- 
ity classes. Nevertheless they are related. The hierarchical renormalization 
group is an ultralocal approximation to a lattice block spin renormalization 
group. See |jlO|, [1^ and references therein. Furthermore, the hierarchical 
model is presumably a zeroth approximation to the full model in a lattice 
derivative expansion. Hierarchical results are to our experience reasonable 
approximations. However concerning their predictions about full models, 
high precision cannot be expected. For instance, u comes out as 0.649 at 
the N = 1 Ising fixed point, as compared to 0.63 from ^-expansion and 
Monte-Carlo methods. Hierarchical models are known to be extremely valu- 
able in gaining qualitative information about the model under investigation 
and to prepare the ground for further renormalization group studies, where 
non-local interactions are included. 

The hierarchical renormalization group for models with an A^-component 
real scalar field is a theory of the non-linear integral transformation 

Ra,p,,iZ){^) = Jdft,{C)Z{P^ + Cr, (3) 

where a, P, and 7 are real parameters, and where dfi^{Q is the Gaussian 
measure on with mean zero and covariance 7 (times the unit matrix). 
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Recall that its Fourier transform is 

d/i,(C)e'^^- = e-^^". (4) 



We consider the renormalization group in the so called high temperature 
picture. See [jl3[ and references therein. Our parameter values in are 



2+D 



a = 2, (3 = 2"^, 7 = 1. (5) 

(You may take this as the definition of the high temperature picture.) Here 
D is the dimension of the model. We restrict our attention to the most 
interesting case, when D = 3. 

The transformation is the composition of three steps: taking the 
square, Gaussian convolution, and rescaling. The numerical computation of 
(H) will be decomposed into these steps. 



3 Algebraic Formulation 

We use the techniques described in ||13|, P). Consider first the 0(A^)-invariant 
case. A general even 0(A^)-invariant Boltzmann factor can be written as a 
sum 

oo 

= ^(02)" Z„, (6) 

?i=0 

and is parametrized by real coefficients Z„. The transformation (|]) becomes 
the following non- linear transformation in terms of the coordinates 

oo m 

R{Z)^ = i?^ ^G„,™(iV,7) (7) 

m=n 1=0 

The coefficients Gn,m{N,'y) are defined as follows. (We call them structure 
coefficients.) Let Pm{4') = {(f^)^- The Gaussian convolution of it is a linear 
combination 

/m 
d/i^(C) + = Y. Gn,m{N,l). (8) 

n=0 
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with coefficients (essentially the number of contractions) 



m—n 



1=1 ' 

It is convenient to define Gn,m{N, 7) = for n > m. 

Analytical (and numerical) experience suggests a different normalization, 
namely Zn = Xn / a/ {2n)\ with a suitable constant p. For notational 
simplicity, we prefer to display the normalization (^). 

For practical computations, one has to truncate the transformation (^ 
to a finite number K of non-zero coefficients. For n < K, the resulting 
transformation is 

K m 

RK{Z)n = G„,„(iV,7) Y,^l^rn-l. (10) 

m=n 1=0 

In the high temperature picture, this truncation scheme is known to converge 



as ^ oo. See |]rT| for a detailed analysis of the one component case. 

Renormalization group fixed points are approximated as stationary fiows 
of the truncated transformation 

RK{Z*)n = K- (11) 

The most interesting datum of a renormalization group fixed point is its 
spectrum, from which one learns the behavior of the linearized fiow in its 
vicinity. The linearized fiows around a fixed point is given by 

K m 

Di?^(F)„ = 2/32«^ G'„,„(iV,7) Y,ZtY„,^i. (12) 

m=n i=0 

The spectrum of Z* is the set of eigenvalues A*^*^ of (|12D. The eigenvalues 
again are directly related to the critical exponents cr'-*-'. See [0. In our 
model, the relation is A*-*-* = a"^^^ . (A peculiarity of hierarchical models is 
that the critical exponents are L-dependent.) The eigenvalue problem reads 

Di?^(r«)„ = A« (13) 

The linearized renormalization group transformation (|^) can be brought to 
a manifestly symmetric form in the no-truncation limit, and is thus diago- 
nalizable. In practice, (|T3|) turns out to be a very reliable way to determine 
the spectrum of renormalization group fixed points. 
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4 Numerical Results: 0(N) Fixed Point 



We have searched for fixed points of the system of algebraic equations ([lOD 
with a Newton algorithm. The program was written in C++ using the data 
type long double representation for real numbers. In order to check roundoff 
errors, we compared the long double with the (simple) double representa- 
tion and found no significant deviations. 
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Table 1: Spectrum of 0(2)3 fixed point 

Table |l] shows the effect of truncation on the first few eigenvalues of the 
non-trivial 0{N) fixed point with N = 2 components in D = 3 dimensions. 
(A trivial "volume eigenvalue" A'-''-' = 2 has been omitted.) We see a rapid 
increase of accuracy with the number of couplings. To get the same absolute 
precision, the smaller the eigenvalue of interest is, the bigger the number of 
couplings has to be chosen. This will be important in the following. 

The first eigenvalue A'-^-* belongs to the critical exponent u = l/a^^\ or 
u = log(a)/D log(A). The numerical value of u is thus 

z/ = 0.7082249. (14) 

With twenty couplings, the first eigenvalue comes out with an accuracy of 
more than three digits, with thirty couplings of more than nine digits. 
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Figure 1: Spectrum of 0(2)3 fixed point as function of truncation order K 
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Table 2: Spectrum of 0{N)^ fixed point at = 30 

A computer assisted proof for the one component case was developed by 
rT| . From it, one gets two sided bounds on the critical exponents with arbi- 
trary precision. It confirms the accuracy of our calculation by a comparison 
at = 1. 

Table ^ shows the spectrum of the non-trivial three dimensional 0{N)- 
fixed point as a function of the number of components A^ at truncation order 
K = 30. 

Analytic continuation in A^ is naturally possible for the system of algebraic 
equations (PUf), since the A^-dependence is encoded entirely in the structure 
constants (H). These depend polynomially on A^. To maintain the same 
precision as A^ increases, one has to increase K. Up to A^ = 3, the truncation 
order = 30 suffices. 

We rediscover among other things the well known result that the theory 
becomes trivial at A^ = —2, i. e. , z/ = 0.5. 
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Figure 2: Spectrum of 0{N)^ fixed point at = 30 



5 Cubic Invariance 



Cubic symmetry is the following. Consider an A^-component model. (The 
local spin takes values in M^.) The cubic symmetry is the finite group of 
transformations, consisting of rotations and reflections, which leaves invariant 
the cube [-1, 1]^. 

We will restrict our attention to even Boltzmann factors, with the prop- 
erty that 

Z((-l)^Vi, • • • , i-iy^'M = Z{cl>,, ...An). (15) 

Such Boltzmann factors are functions of . . . , 0^. To be cubic invariant, 
they have to be in addition symmetric functions of (pf, . . . ,<!)%. An even 
cubic invariant Boltzmann factor is thus a function (pTSj) , which in addition 
satisfies 

Z(0^(i),... ,(/)^(7v)) = Z((^i,... (16) 

for all permutations vr G ©at. We begin with the simplest case, the study of 
cubic perturbations of the 0(A^)-invariant fixed point. 

5.1 Cubic Perturbations of the 0{N) Fixed Point 

We can extend the polynomial basis Pn{4') of 0(A^)-invariants as follows. Let 
n = {uo, Tic) G N X {0, 1} and define 

These symmetric functions do not close under multiplication. They do how- 
ever close under Gaussian convolution. They satisfy 

I d/i(C) PnAi^ + = Y. Prii^) Gn,rn{N, 7), (18) 

with structure coefficients 

G'(n,o),(m,i)(A^,7) = iV 6-0,2(1,7) G'„,m(A^ + 4, 7) 

+61,2(1,7) G'„_i,^(iV + 6, 7), (19) 

G'(„,i),(„,i)(iV,7) = 62,2(1,7) G'„,^(A^ + 8, 7), (20) 
G'(„,o),Ko)(iV,7) = G„,^(Ar,7), (21) 
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where Gn,m{N,'y) is as in (^). The hnearized renormahzation group at the 
0{N) fixed point with cubic 0'^-perturbations comphcates to 

K 1 nio 
2^2n„+4n. ^ G(„„,„,),(^,,^,) (iV, 7) ^(1,0) • (22) 

mo=no mc=0 lo=0 

The 0(A^)-invariants, defined by ric = 0, form an invariant subspace. The 
cubic invariant eigenvectors generally have non- vanishing 0(A^)-components. 
The cubic eigenvalues will be denoted by k,^^\ We order them according to 
their degree of relevance. 
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Table 3: Cubic spectrum of 0{N)^ fixed point at K = 30 

Table ^ shows the 0(A^)-invariant and the cubic invariant spectrum at 
the 0(iV)-invariant non-trivial fixed point in three dimension as a function 
of the number of components N. The leading cubic eigenvector is - in e- 
expansion a deformation of - a cubic 0"^- vertex. Here we restrict our attention 
to eigenvectors of the ^'^-type times powers of 0^. (The subleading eigenvalue 
K^^) belonging to the cubic ^^-interaction is not displayed here.) 

The largest cubic eigenvector, k^^\ becomes one at a critical value Nc of 
A^. We learn from Table | that 2.20 < A^c < 2.25. A closer look at the 
vicinity of Nc yields the following. Table (^ shows that the value of Nc is 
located inbetween 

2.219435 <Nc< 2.219436. (23) 
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(Numerical errors are negligible.) We learn furthermore the important lesson 
that the A^-dependence of k^^^ is rather weak between the two and three 
component models. 
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1.00001973585855 



Table 4: A^'^ at the 0{N) fixed point 



6 Cubic Invariant Fixed Point 

Besides the O(A^) -invariant fixed point, we find a cubic invariant fixed point. 
Again we restrict our attention to the case of three dimensions. The cu- 
bic perturbations (|17D have to be enlarged to a generating system of cubic 
invariant polynomials. We have investigated several possibilities. 
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6.1 Lifted Representation 

The first possibility uses an over complete system. Let n= (ni, n2, na, . . . ) e 
N X N X N X • ■ ■ . Define 

pdt-) = E^' P") 

\a=l / \a=l / \a=l / 

We represent our fixed point by coordinates Z„. With each collection is 
associated a function 

= Pn{<P)Zn. (25) 

||n||<A' 

To define a suitable truncation, we introduce the norm := ni + 2n2 + 
3?7.3 .... In other words, we truncate the model to a maximal power of fields. 
The summation is restricted to the finite subset of given by ||n|| < K. 

Unfortunately, the functions (p^ are not linearly independent. Moreover, 
the linear dependencies vary with N. Therefore, the representation is 
not unique. As an illustration, the situation for iV = 2 is studied in detail in 
the Appendix. 

As we intend to use the over complete representation, we have to specify 
a lift of the renormalization group. We do this as follows. We have that 

Pn(0) P™(0) = Pn+rni(l))- (26) 

The Gaussian convolution can be written in the form 

f dfX,{0 Prni^ + = Yl PnWGn,rniN,^), (27) 

\\n.\\<K 

where the coefficients are the following. We compute 

E 0^^^^'^^ = E ^n,- (28) 
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with structure coefficients given by 



Qn^rn = ^ah Ula rUb 



3 b : < 



Qn.r 



4a ma {ma 



Qn,rn = 2a (2a - 1) ma 



3 a ^ 1 



3 a ^ 1 



(ma -l=na 
mb-l = nb 

ma+b~l + 1 = Ua+b-l 

y c a,b : mc = Tie 
ma-2 = na 

m2a^l + 1 = n2a-l 

\l h ^ a : mb = rib 
ma~l =na 

ma-l + 1 = Ua-l 

\f b ^ a : mb = Ub 



(29) 



gn,rn = 2{2{mi-l) + N)mi 



mi — 1 = 
V a 7^ 1 



ni 
rria 



(30) 



(31) 



(32) 



We then compute the Gaussian integral as the matrix exponential thereof. 
The result is 



Gn,rn{N, 7) = exp Q 



(33) 



The matrix Q becomes upper triangular, when the couplings are sorted ac- 
cording to their total power of fields ||n|| . (When the renormalization group is 
truncated to a finite power of fields, the Gaussian convolution alone does not 
generate higher powers.) Furthermore, the matrix Q becomes nilpotent, and 
the matrix exponential becomes a finite sum. Our first system of equations 
is 



2||n|| 



G„,^(iV,7) 5^ ^/^i . 

||m|j<A' l+k=m 



(34) 



This system of equations defines a lifting of the renormalization group. Every 
fixed point of ( P^ ) becomes through ( pSj ) a fixed point of the original hier- 
archical renormalization group. Their spectra of eigenvalues coincide (in the 
no-truncation limit), but some eigenvalues become degenerate. The eigen- 
value problem is completely analogous to the 0(A^)-invariant case. We do 
not write down the equations here. 
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14742 


8 
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18 


1597 
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18460 


9 


97 


19 
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29 


23025 



Table 5: Number of couplings in the lifted system for given truncation K. 

We can use (|3^) to determine the cubic fixed point for any value of A^. The 
reason is again that N enters polynomially in the structure coefficients (P^. 
The drawback of this representation is that the number of couplings increases 
very fast with the order of truncation. Table (|6.1|) shows this number of 
couplings with ||r2|| < K for K = 1 . . . 29. 

We looked at the system up to = 18. In this model one has to compute 
1597 couplings. To get an idea of the achieved precision, we did one run 
with i^' = 19, which means 2087 couplings. (Has anyone ever computed a 
renormalization group with more unknowns?) 

As in the 0{N) case we recognize one eigenvalue to become marginal at 
a critical number of components Nc. Again the value of Nc is suggested to 
be 

2.2 <N,< 2.25 (35) 

Table (|^) combined with the investigation in polar coordinates for the 
A^ = 2-model, to be described below, suggests that the precision of /i2 at 
= 18 is three digits. 

The calculation with K = 19 was only performed for the case N = 2 
and D = 3, as it required about 22 hours CPU time on an IBM RS6000 
workstation, K = 18 required around 10 hours on the same system. 

As the value fi2 decreases with increasing K, we expect the critical value 
Nc to be given an upper bound by the value Nc < 2.25. 
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Table 6: Cubic spectrum at K — 18. 
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Table 7: Cubic spectrum as function of K 
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To find a cubic fixed point by Newton iteration, one has start from a 
sufficiently close initial guess. We refrain from presenting a lengthy table 
of fixed point couplings. Instead, we present a sufficiently accurate initial 
guess from which anyone can reconstruct the cubic fixed point. (Later we 
will discuss the fixed point couplings in the special case when N — 2). Our 
search strategy was the following. We started at two ends of the interval of 
interest, at = 2 and a,t N — 3. First we looked for the 0{N) fixed point, 
using the start values 

^(0,0,0,0,-) = 0-4 

^(1,0,0,0,...) = 0-06 

^(2,0,0,0,...) = 0-001 
Z{n,o,o,o,...) = Vn > 2. 

Within our scheme (no rescaling on the Z„) this initial guess turned out to 
be sufficient to converge to the the 0{N) fixed point at both values of N. 
To find the cubic fixed point, we took the 0{N) fixed point and changed the 
following coefficients: 

^(0,0,0,0,-) = 0-64 
^(0,1,0,0,...) =0.007 
^(0,0,1,0,...) = 0.00011 



for the case N — 2 and 

^(0,0,0,0,-) = 0.45 
^(0,1,0,0,...) = 0.007 
^(0,0,1,0,...) = 0.00011 



for = 3. Having found the cubic fixed point at these two values of we 
stepped towards the critical Nc from both sides, using the last fixed point as 
starting vector for the next iteration. 
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6.2 AT = 2 Using Polar Coordinates 

For fixed A^, it seems natural to introduce spherical coordinates. For N = 2, 
we have that 

Z{(f)) = Z{r, (^) = ^ Zmir) cos{rrup) . (36) 

m>0 

The cubic symmetry can be implemented by requiring that the coefficients 
Zm{f^) be zero for all m which are not integer multiples of 4. The result Z' 
of a hierarchical RG transformation applied to Z can be written in the form 

Z\<p) = ^-¥'<i>' J d\e-'^^" e^<^^ Z\C) ■ (37) 

We introduce polar coordinates for ^ and through 

^ = r(cos V5, sin v?) , = _R(cos 6*, sin 6*) , (38) 
and expand the square of Z 

Z\0 = Y.{Z'),^{r)cos{mO. (39) 

m>0 

The {Z'^)„i are related to the Z^ through 

i^Z ^fYi 2 ^ ^ ^ ^ Zffi^ ('^m, mi +7712 ^m,\in\—m2\) • (^O) 

mi m2 

The angle integration in eq. (^) can be performed, resulting in 

Z\R,e)=Y,ZM cosimO), (41) 

m 

with 

/■oo 

Z'M = e-^^'""' / t/rr e-^'-' /„,(/3ri?) (Z2)m(r) . (42) 

Here, Im denotes the modified Besselfunction of order m. The next step is 
an expansion of the Zm in a power series of r^, 

Z^ = Y,Zn^y'. (43) 
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Expanding the Bessel function and performing the integration over r yields 
the following relation (for m even) : 



with 



l'-m/2 



A combination of the transformation with the square operation eq. (|40| ) yields 
the structure coefficients of the complete hierarchical RG transformation for 
N = 2. The matrix representing the linearized transformation is then easy 
to compute. 

A Newton solver was used to find the cubic fixed point. It is known 
that for N = 2 the cubical fixed point can be exactly mapped (by a rotation 
of the spin vector with an angle vr/4) to the product of two independent Ising 
{N = 1) fixed points. Our fixed point, plotted in figure ^, is indeed of this 
type. 

We determined the eigenvalues of the linearized transformation at that 
fixed point. Some care was devoted again to study the truncation effects. 
Table |^ shows the six leading eigenvalues at the cubical fixed points for three 
different truncation orders. To the given precision, it is sufficient to include 
K = 32 powers of in the ansatz. Furthermore, it is sufficient to include 
"angular momenta" m with m < rrimax = 4. Going to mmax = 5 did not 
change the results for the exponents at all. From the factorization of the 
cubical fixed point into two Ising fixed points it follows that the spectrum 
can be built from the Ising spectrum through 

Aj,cub = 2 -^i.Isi '^k,M ■ (46) 



The leading Ising eigenvalue is 1.4271725 |T3[. Its square divided by two is 
1.0184107, in nice agreement with the result for the subleading eigenvalue 
of the cubic fixed point. The whole Ising spectrum is part of the cubic 
spectrum as follows when one of the Ising eigenvalues is taken equal to the 
volume eigenvalue two. 

The comparison reveals that the spectrum in the polar coordinate repre- 
sentation and the spectrum in the overdetermined representation agree. Ob- 
viously, K = 18 does not bring about a very high degree of accuracy for the 
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k 


K = 30 


= 26 


K = 20 


K = 19 


K = 18 




1 


1.42717 


1.42718 


1.42867 


1.43034 


1.43377 


1.4337 


2 


1.01841 


1.01840 


1.01831 


1.01844 


1.01889 


1.0188 


3 


0.859411 


0.859393 


0.857925 


0.856580 


0.854073 


0.8540 


4 


0.613255 


0.613058 


0.605915 


0.602012 


0.596873 


0.5968 


5 


0.479634 


0.479563 


0.476318 


0.473871 


0.469780 


0.4887 


6 


0.369252 


0.368630 


0.360018 


0.357220 


0.353718 


0.4697 



Table 8: The first 6 eigenvalues of the spectrum at the cubical fixed point, 
N = 2, for three different truncation orders, mmax = 4 in all cases. For 
comparison the last column shows the results from the calculation with the 
over determined basis, truncation order K = 18, cf. table ( |0|) . 



cubic invariant fixed point in the overdetermined representation. However, 
larger values of K would have required too much computer time. Further- 
more, we learn that the overdetermined representation generates new spuri- 
ous eigenvalues. They should converge to the true ones in the no-truncation 
limit. The first such spurious eigenvalue is however a sub leading (or irrele- 
vant) one. 

It might be interesting to have a look at the fixed point couplings them- 
selves. They are presented for {K, mmax) = (30, 4) in table |[ 

7 Conclusions 

In the framework of the hierarchical renormalization group, we have studied 
the stability of both the 0{N) symmetric and the cubical fixed point for 
D = 3, in the range between N = 2 and = 3. No problems arised 
when investigating the stability of the 0{N) symmetric fixed point. For the 
cubical fixed point, however, the extension of a suitable basis to non- integer 
turns out to be nontrivial. We solved this by using an over-complete set of 
functions. The big number of couplings that had to be used in this approach 
required quite an effort to solve the fixed point equations. Furthermore, it 
might be considered as a problem that the continuation from integer to 
the real domain is by no means unique. The fact, however, that we find 
consistent values of Nc both at the 0{N) and the cubical fixed point seems 
to indicate that the chosen basis is a natural one. One further comment 
is in order: While in the full model the value of A^^^ is very close to 3, the 
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T 

/ 




Zxi 


ry 




ry 





.5668E+00 


.OOOOE+01 


.OOOOE-l-01 


.OOOOE-l-01 


.OOOOE-l-01 


1 


.6404E-01 


.OOOOE+01 


.OOOOE+01 


.OOOOE+01 


.OOOOE+01 


2 


.3162E-02 


.1518E- 


-03 


.OOOOE-l-01 


.OOOOE-l-01 


.OOOOE-l-01 


3 


.9355E-04 


.1190E- 


-04 


.OOOOE-l-01 


.OOOOE-l-01 


.OOOOE-l-01 


4 


.1902E-05 


.4309E- 


-06 


.3827E-08 


r\r\r\r\'\ — \ i r\-\ 

.OOOOE-l-01 


.OOOOE-l-01 


5 


.2877E-07 


.9721E- 


-08 


.2351E-09 


.OOOOE+01 


.OOOOE+01 


6 


.3409E-09 


.1550E 


09 


.6851E 11 


.3149E-13 


.OOOOE-l-01 


7 


.3280E-11 


.1877E- 


-11 


.1269E-12 


.1624E-14 


.OOOOE-l-01 


8 


.2633E-13 


.1807E- 


-13 


.1689E-14 


.4026E-16 


.1135E-18 


9 


.1799E-15 


.1430E- 


-15 


.1728E-16 


.6414E-18 


.5103E-20 


10 


.1064E-17 


.9518E 


18 


.1420E 18 


.7407E-20 


.llllE 21 


11 


.5515E-20 


.5436E- 


-20 


.9673E-21 


.6629E-22 


.1567E-23 


12 


.2532E-22 


.2703E- 


-22 


.5586E-23 


.4798E-24 


.1610E-25 


13 


.1039E-24 


.1185E- 


-24 


.2783E-25 


.2895E-26 


.1290E-27 


14 


.3836E-27 


.4624E- 


-27 


.1214E-27 


.1488E-28 


.8391E-30 


15 


.1283E-29 


.1620E- 


-29 


.4688E-30 


.6628E-31 


.4568E-32 


16 


.3911E-32 


.5135E- 


-32 


.1618E-32 


.2593E-33 


.2125E-34 


17 


.1090E-34 


.1479E- 


-34 


.5024E-35 


.9000E-36 


.8583E-37 


18 


.2786E-37 


.3890E- 


-37 


-I A -t C\~t — 1 O ^ 

.1412E-37 


.2791E-38 


.3042E-39 




/:?rQQT7 /in 

.OOooiIj— 4U 


.9346E- 


-40 


Qcoi^T? /in 
.ooyoiij— 4U 






20 


.1401E-42 


.2046E- 


-42 


.8292E-43 


.1934E-43 


.2633E-44 


21 


.2728E-45 


.4055E-45 


.1719E-45 


.4294E-46 


.6415E-47 


22 


.4759E-48 


.7182E-48 


.3167E-48 


.8407E-49 


.1363E-49 


23 


.7311E-51 


.1117E- 


-50 


.5100E-51 


.1428E-51 


.2490E-52 


24 


.9663E-54 


.1493E- 


-53 


.7021E-54 


.2061E-54 


.3832E-55 


25 


.1068E-56 


.1666E- 


-56 


.8042E-57 


.2461E-57 


.4847E-58 


26 


.9536E-60 


.1499E- 


-59 


.7407E-60 


.2353E-60 


.4879E-61 


27 


.6560E-63 


.1039E^ 


-62 


.5239E-63 


.1721E-63 


.3738E-64 


28 


.3246E-66 


.5174E- 


-66 


.2658E-66 


.8996E-67 


.2039E-67 


29 


.1023E-69 


.1641E- 


-69 


.8570E-70 


.2982E-70 


.7026E-71 


30 


.1539E-73 


.2479E- 


-73 


.1315E-73 


.4694E-74 


.1146E-74 



Table 9: Coupling constants of the cubical fixed point, for K = 30 



20 



hierarchical 0{N) fixed point in three dimensions becomes unstable with 
respect to cubic perturbations already at Ai'c = 2.219. It has been observed 
in other contexts that the dependence of certain quantities on the number 
of spin components is shifted towards smaller values of N in the hierarchical 
approximation. For instance, the two dimensional hierarchical non-linear a- 
model is asymptotically free for > 1 as opposed to > 2 in the full 
setting |l^. This is due to the absence of wave function renormalization 
in the hierarchical approximation. 

Supplementing the present work by a high order e-expansion (which is 
certainly feasible) would be interesting and useful. 
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Appendix: A Basis for the AT = 2 Case 

The polynomials (^) do not form a basis for cubic invariant functions. There 
are linear dependencies. The simplest relation among these is given by 

-P(0,2,0,0) = 2P(o,0,0,l) ~ 2P(2, 1,0,0) + -P(4,0,0,0)! (47) 



where the indices are to be understood in the sense of ( p^ with all rii for 
i > 4 are equal to zero. 

For the case N = 2 a basis for polynomials of type ( PBj ) is given by 



\a=l J a=l 

The relation with ( pij ) is that (^) equals Pm with mi = Uo, m^nc = 1; cind 
all others are zero. 

This basis has the following complications. A product of two such poly- 
nomials is no more the polynomial to the index given by the sum of the 
indices of the factors. In general, a product of two basis elements is a linear 
combination 

PM{mm.,^.)i<P) = E qt£) '^"^°''"^^^(^.,^.)('^)- (49) 

kQ=no+mo 
kc=nc+mc 

We do not have closed expressions for the coefficients Cjj^°'J!^^'^™'°'"^''\ We 
tabulated those we needed in our programs. The table was generated using 
computer algebra. 

The Gaussian integration of polynomials of type (^81) is given by 

d/i,(c) i^K,n.)(^ + c) = E E 4:t;(^^,7)i^(....)(^) + 

ko—O kc—0 f f ni\ 

no n. (50) 

EEc::.^(^'^)^(^.^.-^^)(^) 



ko=0 kc=0 

with 



0";(^' ^) = G^2n.,2..(l, 7)G„„,,„(4n, + AK + iV, 7), (51) 
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where Gn^m{N, 7) is as in (j^). The formula ( ^0] ) does not yet provide a formula 
for the Gaussian integration in the basis (^Sj), as the second term is not 
yet expanded in this basis. Its expansion involves again rather complicated 
structure coefficients which we determined by means of computer algebra. 

We used this representation to check our results with the overdetermined 
representation. This polynomial basis turned out to a less appropriate rep- 
resentation than the following one with polar coordinates. Both suffer the 
drawback that they apply to = 2 only, and cannot be continued to any 
A^. (They can be continued but do not coincide with the desired models at 
integer N.) 
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Figure 3: Visualization of Table 




Figure 4: Cubical fixed point Boltzmann factor Z* for N — 2, 
as function of the two spin components 4>i and 02- 
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